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The strong non-linearity plays a significant role in physics, particularly, in designing of novel 
quantum sources of light and matter as well as in quantum chemistry or quantum biology. In 
simple systems, the photon-photon interaction can be determined analytically. However, it becomes 
challenging to obtain it for more complex systems. Therefore, we show here how to measure strong 
non-linearities via allowing the sample to interact with a weakly pumped quantized leaking optical 
mode. We found that the detected mean-photon number versus pump-field frequency shows several 
peaks. Interestingly, the interval between neighbour peaks equals the photon-photon interaction 
potential. Furthermore, the system exhibits sub-Poissonian photon statistics, entanglement and 
photon switching with less than one photon. Finally, we connect our study with existing related 
experiments. 



PACS numbers: 42.65.Pc, 03.67.Bg, 42.50.Ar, 42.50.Lc 



The non-linear interactions are perhaps the most in- 
vestigated fundamental problems in modern physics and 
related subjects. An enormous amount of novel theoret- 
ical and experimental work has been carried out in this 
field. Significant results re gard ing this issue were ob- 
tained in quantum optics [l-12|, ultracold atomic gases 
U 12 1) quantum electronics [I8I - I27I ] or quantum bio- 
chemistry [28j]. Particularly, strongly interacting single 
photons were discussed in [2j while quantum entangle- 
ment of ultraslow single photons and quantum switch- 
ing were investigated in Ref. 0], respectively. Enhanced 
Kerr non-linearity was reported in (3-01 • Large Kerr 
non-linearity has possible applications in quantum non- 
demolition measurements and for quantum logic gates or 
for optical squeezing and studies of non-locality. The ef- 
fect of Kerr non-linearity on the slow light propagation 
was discussed in Q whereas turning light into a liquid via 
atomic coherence effects was studied in A procedure 
by which strong high-order non-linearities can be syn- 
thesized out of low-order non-linearities was proposed in 
[loj while the quantum dissipative chaos in the statistics 
of excitation numbers was investigated in [ll| . Entangle- 
ment via the Kerr non-linearity in an optical fiber system 
was experimentaly demonstrated in [l2j | . 

A number of experiments were performed emphasiz- 
ing intresting non-linear phenomena in cold atomic sam- 
ples [13Hl6j and trapped single-atom systems in optical 
cavities [171 ]. However, atomic systems are known to ex- 
hibit not so large non-linear couplings compared with 
the corresponding decay rates. Hence, a significant effort 
was devoted towards finding systems showing ultrahigh 
non-linearities. Remarkable, giant Kerr non-linearities 
were found in quantum nanosystems flj| and quantum 
circuit systems [l9j]. The non-liniarities are further re- 
sponsible for a number of fascinating phenomena. Bosc- 
Einstein condensation, bistability and electromagnetic- 
field-induced transparency in high-density exciton sys- 
tems were shown to occur in (20 . l2lj . while polariton 



quantum blockade in a photonic dot was investigated in 
Ref. 22|. Signatures for a classical to quantum tran- 



sition of a driven non-linear nanomechanical resonator 
and photon-number squeezing in circuit quantum electro- 
dynamics were studied in 23[ and 24j, respectively. A 
proposal for detecting singlc-phonon transitions in a sin- 
gle nanoelectromechanical system was discussed as well 
25] . Other non-linear effects involving macroscopic sys- 
tems, superconducting qubits or biochemical samples 
were reported in [26l - l28j further advancing the research 
in these areas. 

In many of the above mentioned effects, at least some 
knowledge on the non-linear interactions is required. A 
number of simpler samples allow the non-linearity to be 
determined even analytically. However, it becomes diffi- 
cult to obtain it for more sophisticated systems. There- 
fore, here, we show how to measure the third-order non- 
linear susceptibility which is responsible for the 
boson-boson interaction potentials, for instance. We 
model the system as a pumped non-linear quantum oscil- 
lator. The examined sample additionally interacts with 
a quantized lossy field mode. The applied field is weak 
such that the non-linearity is not affected by the exter- 
nal driving. We found that in the long-time limit, the 
detected photons show several peaks as function of exter- 
nal field detuning when the photon-photon non-linearity 
is larger than the quantized mode damping. Remarkably, 
the frequency intervals between the neighbour peaks are 
equal with the photon-photon interaction potential al- 
lowing us to determine it. As the strength of the pump- 
ing field is increased, the observed mean photon number 
shows strong asymmetrical behaviors. At lower pumping 
intensities switching with less than one photon occurs. 
Moreover, the switching effect improves for even larger 
non-linear interactions. The photon statistics can be sub- 
Poissonian. These results would disappear for vanishing 
boson-boson interactions. Note that for particular pa- 
rameters one can create an entangled single-photon state 
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|*±) = (|0) ± |1)) /a/2 with a high fidelity. Here, |0) and 
|1) are the number state basis with zero and one pho- 
ton, respectively. Therefore, the study is of significance 
for applications ranging from optical communication to 
quantum processing of information. 

We proceed by briefly introducing the main steps of 
the analytical formalism involved and then rigorously de- 
scribing the obtained results. The Hamiltonian charac- 
terizing the interaction of a non-linear oscillator possess- 
ing the frequency lo with a coherent source of frequency 
ojl, in a frame rotating at ujl, is: 



H = -Aa'a T aa^a 2 + e{a) + a) 



(1) 



where A = lul — lo and a > signifies the photon-photon 
non-linearity (a oc Refx*- 3 -*]), while e is the amplitude of 
the coherent driving, and a are the creation and the 
annihilation operator of the quantum oscillator, respec- 
tively, and obeying the standard bosonic commutation 
relations, i.e., [a,a^] = 1, and [a, a] = [a^,a^] = 0. The 
=F sign in Eq. (QJ accounts for an attractive or repul- 
sive boson-boson interaction, respectively. Actually, the 
Hamiltonian (|TJ) is abundantly investigated and charac- 
terizes a wide range of processes. A formidable approach 
to analyze the driven damped non-linear oscillator is via 
the master equation. Hence, the system is described by 
the reduced density operator, which in the interaction 
picture and under the usually applicable Born-Markov 
and rotating-wave approximations satisfies the master 
equation [iillicjl]: 



P = -i[H, p] - n(\a\ap] + [pa\a]), 



(2) 



where k is the decay rate of the quantized optical oscil- 
lator mode and the overdot denotes differentiation with 
respect to time. 

In order to investigate Eq. ([2J) in the long-time limit 
and for weak excitations in more details, we first shall 
apply the Holstein-Primakoff transformations [3l|. and 
then follow the solving procedure developed in [32j for 
two-level particles. We define a state |g) denoting N — q 
spins pointing up and q spins down. The raising and low- 
ering spin operators act on this state as follows: S + \q) = 
y/q{2s -q + l)\q - 1) and S~ \q) = ^{q + l)(2s - q)\q + 
1), where 2s gives the number of excitations (spins) 
in the system, and < q < 2s. The spin operators 
satisfy the usual commutation relations of su(2) alge- 
bra, i.e., [S^S*] = ±S ± and [S+,S-] = 2S Z , where 
S z \q) = (s — q)\q). Taking into account these properties, 
one can deduce that [31|: 
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Substituting ([3]) in the steady-state form of Eq. @, ob- 
serving that S~ S + = 2sa'a — a) a 2 , and assuming lower 
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FIG. 1: (color online) The mean photon number (n) as well 
as the normalized second-order correlation function g^ 2 \0) 
as function of scaled detuning A/a. Here (a,c) fl/a — 
0.06,0.04,0.006 while (b,d) ft/a = 0.3,0.2,0.1 for the solid, 
the dashed and the dotted curves, respectively. Other param- 
eters are: j/a — 10~ 3 and 2s = 50. 



excitations, i.e. a* a/2s <C 1, one arrives at the following 
long-time solution for the normalized diagonal elements 
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(q\p s \q) of the density matrix p s : 
(q + m)\(2s 
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Here, the normalization constant Z is given by the ex- 
pression: 



2 s 



n 2 



m=0 



(2s + m+l)!(m!) 
(2s-to)!(2to + 1)!' 



(5) 



- 2m \T(l + m + i4>*) /[m\T{l + i4>*)]\ 2 and 
±(2sa ± A)/ (7 ± ia), while Q = 



where C mm = \a 
a = 0/(7 ± ia), (j) 

e/V2s and 7 = k/2s. With the help of the steady-state 
solution (j4]), one can calculate the oscillator's variables of 
interest. In particular, the mean number of excitations 
in the system (n) = (a^a) as well as its second-order 
correlation function G^ 2 '(0) = 
from the folowing expressions: 



t 2 

7 I 



a ) can be obtained 



2s 
q=0 



2.s 



C7( 2 )(0) = ^ O ( O -l)P 9 . (6) 

q=0 



Because we considered the weak excitation limit, the so- 
lution Q is valid as long as f2 -C a. This means that 
the amplitude of the external driving e should be less or 
of the order of boson-boson potential a. Thus, our ap- 
proach is well suitable to describe low-photon processes 
with mean-photon numbers of the order of few photons. 

We further focus on investigating the properties of the 
driven and damped non- linear oscillator using Eqs. (131111) ■ 
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As a first result, Fig. (p} shows the steady-state depen- 
dence of the mean photon number (n) and the normalized 
second-order correlation function g^(0) = G^ 2 \0)/(n) 2 
versus the external field detuning A/a for the nega- 
tive sign in Eq. (P). There, (a,c) describe the weak 
field limit whereas (b,d) depict pumping with a mod- 
erate coherent source. Interestingly, single-photon light 
with sub-Poissonian photon statistics can be obtained 
(see Fig. [Ha, c) near A/a = —1). This is not trivial since 
even weaker fields lead to super-Poissonian photon statis- 
tics (see Fig. [IJa,c) at A/a < —1). However, the two- 
photon correlator (0) has small values here and larger 
values for g( 2 \0) are due to even smaller denominator, 
i.e., (n) 2 . For stronger driving, super-Poissonian pho- 
ton statistics almost always occurs for A/a < —1.5 and 
Poissonian or sub-Poissonian statistics for A/a > —1.5 
(see Fig.[TJb,d)). Furthermore, the mean photon number 
shows strong asymmetrical behaviors. The asymmetry is 
more pronounced for intenser external fields (compare 
Fig. QJa) and Fig. H{b)). Notably, the multi-peak be- 
haviors of the mean photon number can help us to de- 
termine the non-linearity a. One can observe that the 
scaled frequency interval between the neighbour peaks 
in Fig. [lja,b) equals unity and, thus, the corresponding 
inter-peak frequency will be equal to the non-linearity 
a. Therefore, the boson-boson non-linearity can be ex- 
tracted by measuring the mean photon number against 
the external pumping field detuning. The only condi- 
tion is that the non-linearity a should be at least few 
times larger than the damping of the quantized mode 
k. This is well-satisfied in a wide range of systems fl9| . 
Note that the peaks at A/a = —6 in Fig. ((T}d) and 
Fig. (|T}1) are slightly smaller than those obtained numer- 
ically. The inter-peaks interval is obtained correctly for 
even stronger pumping, however, the magnitude of the 
peaks as well as their positions will be given only aprox- 
imately. The same results occur for the positive sign 
in Eq. ([T]), i.e., for repulsive boson-boson interactions, 
though the corresponding curves in Fig. ([1]) will be mir- 
rored with respect to the 0Y axis. Therefore, for simplic- 
ity, we shall consider further only attractive boson-boson 
interactions. 

An explanation of the multi-peak behaviors in Fig. (pj 
can be found by representing the Eq. ^ via the number 
state basis, i.e. P mn — ("t-IPsI"-)- Then, for the negative 
sign in Eq. ([TJ, one arrives at the following terms next 
to P m n, that is, P mn [i(m — n){A + a(m + n — 1)} — 
K,(m + n)]. One can observe here that resonances occur 
at A + a(m + n — 1) = 0, when a > k. The off-diagonal 
elements induced by the driving field in the steady-state 
play a crucial role here. In their absence the mean photon 
number, for instance, would be insensitive on photon- 
photon interactions. In a recent work |33j |. additional 
resonances were found which may help to extract the non- 
linear parameters, though in a more complicated setup. 

To further describe our system, we plot in Fig. @ 




2.0 F 



(b) 



1.0 -/■ 
I 



0.00 0.05 0.10 0.15 0.20 0.25 

Q/a 



0.5 A 



0.0 L 



0.00 0.05 0.10 0.15 0.20 0.25 

Q/a 



FIG. 2: (color online) The mean photon number (n) versus 
the intensity of the pumping field Q/a. Solid line is for A/a = 
—2, dashed curve corresponds to A/a = —1 while the dotted 
line to A/a = 0. Here, (a) -y/a = 1CT 3 and (b) j/a = 1(T 4 
with 2s = 50. 



the dependence of the non-linear oscillator's mean pho- 
ton number as function of the applied field intensity. 
The detunings are adjusted to values giving the peaks 
in Fig. ((T^). Switching with less than one photon is ob- 
served here (see the dotted line in Fig. [2J. The mean pho- 
ton number abruptly increases as the pumping strength 
varies only a little. For instance, the transition among 
states with (n) = to (n) rj 1 occurs while the in- 
tensity of the external pumping field changes between 
Q/a ~ 0.005 and Q/a « 0.02 (see the dashed curve 
in Fig. [2b)- The switching phenomenon considerably 
improves for larger boson-boson interactions (compare 
Fig. [2^ and Fig. [2]d). Positive values of A, away from 
resonance, do not lead to switching effects. These criti- 
cal behaviors take place for repulsive photon-photon in- 
teractions too, but for A > 0. Thus, the appropriately 
prepared non-linear sample can be presented as a single- 
photon optical switching devise. Finally, for comparison, 
we have obtained the oscillator's mean photon number 
and its second-order correlation function without tak- 
ing into account of the boson-boson interaction potential, 
i.e., we set a = 0: 



G^(0) 



(k 2 + A 2 ) 2 



(7) 



In this case, however, the rich variety of effects shown in 
Fig. ([T]) and Fig. (|2|) vanishes and the photon statistics 
allways shows Poissonian photon distribution. 
In the experiment described in Ref. 



15j, few-photon 



switching, Kerr non-linearity and dispersive optical bista- 
bility of a Fabry-Perot optical resonator due to the dis- 
placement of ultracold atoms trapped within were re- 
ported. There, the photon asymmetry similar to the one 
shown in Fig. QIa,b) was observed as well (see also Fig. 2 
in (l5j] ^ - Though experimentally, they got multiple peaks 
of the intracavity mean photon number their nature was 
not explained. Based on our results, we can conjecture 
that the interval between the main peaks around the cen- 
tral one of the intracavity mean photon number obtained 
in [15[ may give the induced boson-boson non-linearity. 
Optical switching with single photons was observed as 
well [341. 
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cess. Sub-Poissonian photon statistics occurs for particu- 
lar external controlable parameters. The reported results 
can be tested with existing experiments. 
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FIG. 3: (color online) The mean photon number (n) as well 
as the fidelities 4>± against the pumping field detuning A/a. 
The solid line depicts the mean photon number (n) while 
the dashed and the dotted curves correspond to <E>_ and $+, 
respectively. Here, Q/a = 0.06, j/a = 10 -3 and 2s = 50. 



In order to look at more details in the few-photon pro- 
cesses emphasized here, in Fig. ([3]), we plot the fidelities 
$± = (*±|p s | v E'±) of entangled states 
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|* ± ) = - 7 =(|0)±|1)), 



(8) 



where |0) and |1) are the number state basis with zero 
and one photon, respectively. There, the dashed curve 
corresponds to $_ while the dotted one to The solid 
line shows the non-linear oscillator's mean photon num- 
ber (n) . The parameters are the same as for the solid line 
in Fig. (pji). One can observe here that at A/a = 0.5 the 
entangled state |^ + ) = (|0) + 11))/-^ is created with a 
fidelity $ + w 0.9 and (n) w 0.4, while the photon statis- 
tics is sub-Poissonian (see Fig.QJ). On the other side, at 
A/a = -1.5 the entangled state |*_) = (|0) - |l))/\/2 
occurs with a fidelity $_ « 0.7 and (n) « 0.2, and where 
photon statistics is super- Poissonian (see Fig.[TJ;). Larger 
fidelities mean establishing of quantum coherences in the 
system and, thus, proving the quantum nature of the 
states ©. Note that for detunings approximately within 
-2.1 < A/a < one has $ + + < I denoting the 
existence of higher Fock states, i.e. \n > 1), though with 
smaller probabilities (see Fig. [3]). Therefore, we have 
demonstrated here switching and entanglement with less 
than one photon processes. 

In summary, we have shown how to measure the 
photon-photon interaction potential via weakly pump- 
ing the quantized mode with which the examined system 
interacts. The non-linear effect induces multiple peaks 
in the mean photon number steady-state behaviors with 
inter-peaks intervals equating the required non-linearity. 
For this to occur the boson-boson non-linearity should be 
larger than the quantized mode damping. This require- 
ment is well-fulfilled in a wide range of systems. Due to 
the involved non-linearity, an asymmetry in the first- and 
second-order field correlation functions was observed as 
well. In addition, we have demonstrated switching and 
high-fidelity entanglement with less than one photon pro- 
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